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method at low and intermediate temperature is proved as well as its inadequacy to disclose quantum 
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1 Introduction 

The thermodynamics of most everyday- scale phenomena, e.g phase transitions, often reflects the un- 
derlying existence of a complex microscopic many-body dynamics which is approachable, generally 
speaking, in statistical terms only. Understanding the collective behavior of physical systems compris- 
ing a myriad of interacting or not, classical or quantum particles does constitute an endless challenge, 
over the last years fed by the astonishing progresses in the production of new materials having prop- 
erties tailored on demand [Sl fT^I^ , as well as in the manipulation of mesoscopic quantum systems 
119] . Experimental and technological successes of this kind do indeed open new problems on a variety of 
many-body systems [51ll6| , to be dealt with by making recourse to the ideas and methods of statistical 
physics (see for example [311] )• The correspondent task of the theoreticians attracted by these prob- 
lems is to contrive simple enough but well addressed microscopic models as well as to develop, within 
the Gibbs framework, new widely applicable mathematical/numerical tools leading to predictions of 
experimental interest. Numerical approaches based on the Renormalization Group (RG) method [SI 
I22j have been successfully employed to produce a deeper insight into the behavior of quantum mat- 
ter. Such applications are even more important in the cases of non-integrable models when an exact 
(even if usually complicated) solution is not available. Nevertheless these methods often require a huge 
computational effort, and might thus be heavily computer-based |17j . 
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Recently a different approach based on a generalized version of the Central Limit Theorem (CLT) 
has been proposed by Hartmann, Mahler and Hess [TU] and, at least for some specific class of problems, 
it turns out to be much simpler in its application than the RG-based treatment. Such an alternative 
method relies on a gaussian approximation for the distribution of physical properties of quantum 
systems and in what follows we will refer to it as the Gaussian Approximation (GA). In this paper 
we critically review such a method, clarifying its limits and highlighting its range of applicability, also 
providing applications to some exemplary physical system. In particular wc show that this method has 
good agreement with the exact thermodynamic behavior only in the limit of high temperature, while 
it fails in reproducing the physics of a wide class of systems at low and intermediate temperatures. 
The existence of limitations for this method has already been noticed in [T^ , where the authors show 
that GA is unable to predict the true behavior of a particular system, namely an Ising chain, at low 
temperature, and suggest that this should be a generic feature of GA. The novelty of our work is that, 
for the first time, it is rigorously shown the existence of systematic drawbacks in the applicability of 
GA. In particular our analysis confirms, elucidates and generalizes the limitations of this approach as 
discussed in [TU] . 

This paper is organized as follows: Section 2 is devoted to a brief review of GA. Mathematical 
techniques allowing for closed expressions for thcrmodynamical quantities in the framework of GA are 
developed in Section 3. The method is applied in Section 4 to obtain results on the thermodynamics 
of spin models, paradigmatic examples of many-body quantum systems. Comments and conclusions 
are drawn in Section 5. 



2 Gaussian Approximation 

Here we briefly review the Gaussian method, originally presented in jlO] and based on a mathematical 
theorem proven in [11] . The notations employed strictly follow [10] . For more mathematical details on 
the fundamental theorem underlying the approximation and its proof wc refer the interested reader to 
[TT] and to references therein. 

Let us consider an A^'-block quantum linear chain with Hamiltonian H = Xl^=i where is 
an operator composed of the block self-Hamiltonian and the interaction between blocks n and + 1, 
thus reading 

I being the identity operator in the single-block Hilbert space. Each block is described by a finite- 
dimensional Hilbert space with dimensionality d, the dimension of the total Hilbert space thus being 

. It may consist e.g. of one or several neighboring spins. Given a total Hamiltonian H , one may 
define the partitioning into blocks in many different ways. It is worth stressing, however, that GA 
applies as long as the dimension of Hilbert space of each block stays finite. 

Let us now define two basis of the total Hilbert space of the system. The first one, whose vectors 
will be written as is composed by eigenstates of H such that H\(f)) = The second basis is 

made up of factorized vectors, written as \a) = I^/j); where each \a^) is a local state of block /.t. 
Finally, we define the two Gaussian parameters 

E-a^{a\H\a), 

Al = {a\H^\a) {a\H\af . 

To calculate free energy and other thermodynamic quantities in Gaussian approximation, one needs 
to know only these two parameters and the ground state energy. Let us for a moment leave apart 
the problem of calculating the ground state energy (which may be solved sometimes exactly and 
otherwise numerically). It is often an easy task to calculate the above two gaussian parameters, as 
will be demonstrated on some examples below. For this reason Gaussian method may seem apparently 
more attractive than other, much more sophisticated methods. 

From the definitions given above, it is possible to build the operator 



(3) 
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with eigenvalues z^. It was shown in with the use of a Central Limit Theorem, that under the 
condition that a constant C > exists such that 



Al > 7VC, (4) 



the following limit strictly holds: 



lim Va{z^ £[zi,Z2\) = I —=dz, (5) 



where Pa(z0 G [21,22]) = E{|0>:zi<^^< 



^2} 



ah 



If one naively applies this result for calculation of density of states rjQ{E) and partition function 
Zo, one readily obtains the following expressions [10], which constitute what we call the Gaussian 
Approximation: 

where Eg is the ground state energy of the system, the summation is extended over the whole factorized 
basis {|a)} and erfc(a;) is the conjugate gaussian error function defined as 

erfc(a;) = —= \ e^^'dw. (8) 

Eq. (O is obtained from eq. ([5]) via 

/.+00 

Zg= / TjG{E)e-^''dE. (9) 
Here /? is the inverse temperature in units of fc^ . 

It should be emphasized that rigorously speaking eq. ([5]) does not imply eqs. (|6l), Q. The reason is that 
the convergence in ([S]) is the convergence of a distribution, which is too weak to justify the summation 
of large (infinite, in the limiting case) number of terms [TO]. Thus eqs. ([7]) must be considered as 
an approximation motivated by a Central Limit Theorem, whose range of validity should be carefully 
established. This latter task is the goal of our work. 

In performing such an analysis, we often will work in the limit of very high number of blocks, 
when the number of terms in the sum in eqs. ([6]), ([7]) becomes very large and goes to infinity in 
the thermodynamic limit. In these situations it is reasonable to exchange summation over {\a)} by 
integration over energy, introducing a density rjfreeiE) defined as the number of states \a) for which 
Ea G [E,E + dE]. With this definition one can then replace summation with integration as 



E 



{E)dE. (10) 



3 GA and exact physical predictions: a comparison 

In this Section we are interested in studying the physical predictions of GA and comparing them with 
the true behavior of thermodynamic quantities. In what follows we try to be as general as possible, in 
order to analyze the how and why GA shows limits of applicability. In the next section we will show 
how these limits apply to some particular system. 
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3.1 Error function at large positive and negative arguments 

To begin with, we highlight a technical issue which is crucial for discussing GA at large N. Let us 
focus on the behavior of complementary error function which enters cq. ([7]). Its argument grows with 
number of blocks as \/N [T^]. Indeed, one expects that 

E^ = Nsg{l + 0{N-^)), 

Ea=Nea{l + 0{N-^)), (11) 
Al=N6l{l + 0{N'^)), 

where Eg, Sa and 5a are quantities defined per block and thus do not depend on N. The above linear 
scaling is a generic feature of a modular Hamiltonian ([T|) . This will be exemplified by means of specific 
models in Sec. |H Also note that linear scaling of is in accordance with the condition (0]). So the 
argument of the complementary error function in eq. ([7|) can be written as 

^ eg--a+p5l ^ (12) 

At large N this is a large number, either positive or negative depending on the sign of (e^ — + j35'^) 
[12| . Complementary error function has a well-known asymptotic behavior at large arguments: 



ertc(a;) ^ \ -x'^ ,- , (13) 



As we will see this different behavior will produce certain discontinuities at the temperature where 
(sg ~ £a + /3(5^) changes its sign. In what follows we will need not only the above asymptotic relations, 
but also an exact inequality [T] which limits erfc(a;) at x > : 

erfc(2;) < e-^'/lv^^;). (14) 



3.2 High temperature limit 

To begin with, we want to compare the high temperature behaviors of both the exact partition function 
Z = Tr(e^'^^) and the partition function obtained from GA. Formally expanding the exponential 
function of H in Z up to the second order in /3 one straightforwardly gets 

Z = Tv(l- I3H +^hA +0{I3^). (15) 



Note now that such a trace is independent on the basis chosen to practically evaluate it. Thus one can 
use the factorized basis {|a)} defined in the previous section. In this way, equation (jisp can be written 

as 



{ia>} 



2 2 (16) 

^ (l - mm + y (a|H|a)2 + ^{{a\H'-\a) {a\H\af)) + 0{p\ 



Equation (|16p does not of course coincide with ([7]). Nevertheless, in the limit of high temperature, 

T > max _ ° , (17) 

{\a)} Ea — Sg 
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arguments of all error functions in ([7]) are positive, and one can use the asymptotic form of erfc(a;) at 
X — > — oo, sec cq. (fT5)) . to obtain 

{|a>} 

= ^ (l - mm + ^m\o^? + ^{{a\H^\a) {a\H\af)) + 0{li% 

which is the same resuh obtained in (|16p . 

It is thus shown that GA and exact resuhs coincide, up to the order in the high temperature 
hmit, since by exploiting the knowledge of partition function one can obtain expressions for every 
thermodynamical quantity of interest. 



3.3 Low temperature limit 

What about low temperatures? In order to avoid unwanted divergences when taking the limit T — >■ 0, 
in the present subsection (but nowhere else!) we define the Hamiltonian ^ in such a way that the 
ground state energy Eg equals zero. Then the exact partition function reads 



Z = l+^e-'3^' > 1. (19) 

i=2 

Here Ei are eigenstates of the total Hamiltonian. 
On the other hand, at low enough temperatures, 

r<min-2-, (20) 

{\a)}ea 

one is able to limit the approximate partition function Zq from above. Indeed at such low temperatures 
arguments of all error functions in eq. ([7]) are negative, and exploiting inequality (|14[) one obtains 

Zg < max , ^0 as T ^ 0. (21) 

{|a>} y^^d^{l36a-ea/6a) 

We thus see how the approximate partition function can not reproduce the exact result when temper- 
ature is low enough. This is a consequence of the fact that the Gaussian Approximation does not hold 
when the system is very close to being in a single eigenstate of H (in our case, in the ground state) 
|10j . In jlOj such a failure in reproducing exact results at T « was noted in connection with the study 
of a Ising chain, and it was argued that GA would in general fail when predicting physical properties 
at zero temperature. 

Up to now we have conducted a very general analysis of GA at high and low temperatures, which 
has revealed that while at high temperatures this approximation works satisfactory, at low tempera- 
tures something goes wrong. In the next sections we show that under specific constraints on gaussian 
parameters Ea and Aa it is possible to give explicit expressions for measurable quantities obtained 
from Zg in the thermodynamic limit in the whole range of temperatures. This analysis will allow us to 
elucidate that, contrarily to what happens in the high temperature limit, GA is unable to reproduce 
physical results in the range of low and intermediate temperatures. 
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3.4 An analytical expression for the free energy under constant Gaussian parameters 

Consider a model in which Gaussian parameters Aa and Ea in some factorized basis do not depend 
on a: 

Aa = VNS > V a. 

Note that the first equation implies the ground state energy to be negative unless the Hamiltonian is 
identically zero. From equation ([7]) it is then possible to obtain the Gaussian Approximation for the 
density of states and the partition function: 

_ " -E^/{2S^N) 

\T S J 



Za ~ — e 2 erfc 
2 



We wish to calculate free energy per particle which is defined as — iV^^Tln Z. In the GA framework 
it can be easily obtained from eq. (j23p with the use of asymptotic expressions (jl3p for error function. 
In the thermodynamic limit the result reads 

la. li,n(-iV-.T,„Z„)J-^0"''-A+^»'«'-<-f' (24) 

y-T\nd-^ forr>-|^. 

Let us compare the low-temperature behavior of the exact free energy / and the approximation Jq . 
First of all note that, at T = 0, /g = Sg as it should. One can easily verify that it is the choice of 
ground state energy as a lower limit of integration in eq. ^ which ensures that fo at zero temperature 
is always equal to the ground state energy per spin. 

What about the entropy at zero temperature? GA gives Sg = ^^|t=o = — lnd+ However, 
according to the third law of thermodynamics, entropy is zero at T = [14|FI This disagreement shows 
how GA is unable to reproduce exact thermodynamics when temperature is low enough. 
Note further that 

— Jg and are continuous at T = — |-, however ^yf- and further derivatives are not. Thus GA 
predicts a fake thermal phase transition. This shows how, even at intermediate temperature, the 
agreement with exact results is very poor. 

— /g at r > — J- coincides with the first two terms of high-temperature expansion of /, as expected 
from (E 



3.5 A generalized formula for the free energy under weaker constraints on Gaussian parameters 

Let us now relax one of the two constraints imposed on Gaussian parameters in the previous subsection 
and let us consider the situation in which Aa = const but Ea is not. 
With the help of (fTO)) we can rewrite eq. ([6]) in an integral form: 

m{E) = J ,/,.ee(^)-^e-(^-^)V(2^^^)d^. (25) 

Here the "free" density of states -qfreeiE), introduced in (jlOp . is itself expressible as a summation 
S{|a)} ~ Ea)- It is often useful to approximate it by some more regular and manageable function. 
In what follows we approximate it as a gaussian function (as often done when dealing with these kind 
of densities [13]), exploiting a second time the Gentral Limit Theorem on which GA is based, thus 

^ The third low can be violated only for systems whose ground state degeneracy grows exponentially with 
the number of blocks. 
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introducing a further source of error since the gaussian is unable to reproduce exactly the behavior 
of r]free{E) at its tails. We call the approximation thus obtained Double Gaussian Approximation 
(GGA), stressing the fact that the CLT is applied twice in two different steps of its derivation. Being 
aware of the further error thus introduced, let us develop the GGA. Let us assume that single-block 
free Hamiltonians ff^ are of equal form for all (generalization to the case of different i/^ is straight- 
forward). Consider a basis {|aji)}, I = 1,2, d in a d-dimentional single-block Hilbcrt space. Without 
loss of generality we may put 

I. (26) 



The statistical standard deviation of a single-block sclf-Hamiltonian relative to a specific basis {|ajj)} 
is then defined as 



1 



d ^ 

I 



From the Central Limit Theorem we thus obtain 



rifreeiE) 



V27riVcr 



(27) 



(28) 



Note that, as already pointed out, such an approximation is not a good one at the tails of the gaussian. 
However exploiting it, and within its limits, it is possible to apply the GGA to obtain 



Vg{E) 



' - d^ 
dE- 



^-{E-Ef/{2S'^N)-E^/{2a^N) 



2T:NaS 



^N^^e^N/{2{S^+a^)) 
2y/2TTN{S^ 



erfc 



£l - 



(29) 



Here E = Ne and ei 



min (a Ji/^|a J 

{|a;,>} 



Note that the case studied in the previous subsection, namely 



Ea = const, corresponds to the limit ct — > 0. In this limit erfc(x) ^ 2 and eq. (|29p actually reduces 
to the expression given in (j23|) . as it should. In the framework of this approximation and using the 
asymptotic expressions for erfc one obtains 



N 



1+^ 



--l3eN 



2^27rAr(J2 ^ ^2) 



1+^ 



de + N 



-PeN 



y/27rN{S^ + cr2) 



-de. 



(30) 



where 0{x) is the Heaviside step function and = max|eg,ei(l -I- f2)|- 

The first term corresponds to the erfc(a;) with positive argument, while the second term corresponds 
to the erfc(a;) with negative argument. There are some cases where the first term can be neglected or 
equals zero (e.g. in Ising model). The second term alone gives rise to another complementary error 
function which closely resembles the one analyzed in the previous subsection in order to obtain ca. (j24p . 
Following the same steps we can then get 



-T\nd- 



2T 



))+e„ forr< 

for T > -*-t2_ 



(31) 



Equation (pij) shows the same behavior as P^ . All comments made for are thus valid also in 
this case. 

It is worth stressing that ([5l]) is valid as long as the first term in ([50]) can be neglected or actually 
equals zero, and within the limits of the additional approximation used to obtain 
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4 Application to specific models 

As an example of application of the Gaussian method to physical systems, we study its predictive 
power for the quantum Ising chain, extensively studied in literature [T51[^. 



4.1 Gaussian Approximation applied to the Ising model 

Quantum Ising model describes interacting qubits embedded in an external transverse magnetic field. 
It is probably the simplest quantum integrable model available. Its Hamiltonian reads 

N N 

Hi, ^-JYI +''^^h (32) 

where u", a = x,y,z is the a-component of the Pauli matrix vector corresponding to site j, J is the 
spin-spin interaction parameter and h is the external magnetic field. In what follows we mainly will 
work in the thermodynamic limit, where different boundary conditions yield the same result. 

Being integrable, this model can be exactly diagonalized via a Jordan- Wigner transformation |15j . 
obtaining the following expressions for ground state energy and for free energy per spin in the ther- 
modynamic limit 



fis 



7r 

- / dpJj"^ + h'^ + 2hJcoap, 
71- J 



TT 



T 

TT 



dp In 



, a/J2 + + 2/iJcosp 

2 cosh — 

T 



(33) 



(34) 





Let us introduce the following notation for the eigenstatcs of cr^ 

^''\y) = \y), = 

and analogously for x and z. 

In this notation |z) = | t), |^) = I i)- 

To apply GA, we take as factorized states the eigenstates of cr'^ , e.g. 

\a.y) = \yyyyyyyyy--y), (35) 

easily obtaining Ea ~ and 5^ ^ + h^. We now may apply formula (p4)) to obtain f^, which reads 



^"-|/|+forr>-4&^; 



f^'- for T < 



where 



f-- =^T\n2- —-^^^——^( I dp^/j^ + + 2h J cos p] --[ dpy' + Ii^ + 2hJ cosp, 
fl;+ = -T\n2- 



2T 



(37) 



As expected, this equation predicts a wrong behavior at low temperature, as can be readily checked in 

Fig.m 
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4.2 Gaussian Approximation for the Ising model in a different factorized basis 

Let us now consider a different factorized basis, namely the eigenbasis of iJ*' ^ hlZ^=i'^j j which is 
the same basis used in flOj to numerically evaluate the partition function from equation ([7]) . With the 
same notation introduced in the previous subsection we set 

\az) = \zzzzzzz...z). (38) 

This gives 

El = —h, 5z = J, cTz = h (39) 

and 

'E'a = h{2k- N), (40) 

where k is the number of up spins in the factorized state ja^) or, which is the same, is the number of z 
in the state \az). Ea in this case actually depends on a. Equation might thus be applied, provided 
its conditions are met. 
One can check numerically that 

e,il + Sl/al)-el^<0 (41) 
for all values of model parameters J and h. Thus free energy is given by eg. pip . However, as can easily 
be checked, this result exactly coincides with one obtained in the previous subsection employing the 
y-basis. The reason is that when we take a y-product basis as done before we get 

Sy = V J2 + /i2 (42) 

and as a consequence 

5l = 5l+al (43) 

This relation, here shown for the Ising chain, indeed holds for a certain class of systems as expressed 
by the following 

Statement: Consider a model with Hamiltonian ([T]) represented a.s H = + 1 . I being an inter-block 
interaction, and two product bases, {|a)} and {|a)}, such that 

1. Aa = const and Aa = const, 

2. Va Ea = 0, 

3. Va {a\I\a) = and H°\h) = El\h). 
Then 

~5^ + ~a^ = 6\ 

The above Statement implies that if conditions 1,2 and 3 are met for a particular system, then 
equation and equation ([5T|) actually give the same result and one is entitled to use the simpler 
one. We might have stated this Statement right after equation (|5T|) in the Section devoted to general 
results. However, for the sake of clarity and to better show its meaning, we decided to formulate such 
a Statement together with a concrete example. 



We plot and obtained both from a GA in y-basis and from a GA in z-basis, on Fig. [TJ 
One can notice a very good agreement between both the approximate free energies and the exact one 
as long as temperature is high enough (T ^ ^^)i expected from the high temperature expansion, 
equation (jl8|) . but such a good agreement breaks down at intermediate and low temperatures. GA 
in y-basis, in particular, gives the worst approximation, while GA in z-basis is closer to the exact 
free energy, while still showing all the drawbacks highlighted in the previous Section. Note that GA 
in y-basis, as shown, actually coincides with GGA for this particular model and thus is affected by 
a greater source of error, while GA in z-basis, as derived directly from ([7]) (and as such, not in the 
thermodynamic limit), is more accurate. The difference between the approximate and exact results is 
even more marked if we consider that all measurable quantities are obtained as derivatives of the free 
energy with respect to temperature. GA thus predicts, even in the thermodynamic limit, a behavior 
which is, at low temperature, very different from the one actually shown by Ising chain. 

Note also that GGA can be applied in y-basis also (with 5y = h), however the dependence of Ea 
over I a) is more complicated in this case. 
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Fig. 1 Free energy per spin plotted against temperature obtained from the exact formula (|33|) (green solid 
line), from GA in y-basis (obtained both from p6|) or from the same expression stemming from H31|) . black 
dotted line) and from GA in z-basis (obtained directly with the use of ((T]) in the case of A'' = 1000, red dashed 
line) for an Ising chain when J = 1 and /i = 2. It is possible to notice, as expected, how the high temperature 
behavior predicted by GA closely resembles the exact one, while at low temperature GA in both bases violates 
the third law of thermodynamics. Nevertheless the agreement with the exact result is better for GA in z-basis. 

4.3 Other models 

In the present subsection we give examples of other models to which the GA and GGA can easily be 
applied (showing, however, all the drawbacks discussed above). 

— cluster Ising chain PO] : 

Hci = ~'JYl,j=i ^J-i^l^J+i + ^ fj-ifj • In this case one can apply both GA according to eq. 

(in X- or z-basis) and GGA according to eq. (|5T|) (in y-basis), and easily check that the three 
conditions of the Statement actually hold. Thus GA and GGA coincide in this case. This fact, and 
further the fact that the application of GA yields exactly the same result as the one obtained for 
the Ising chain in the previous section are not surprising: it has indeed been shown that the Cluster 
Ising model is equivalent to an Ising chain in external magnetic field, and as such also GA can be 
applied equivalently. 

— cluster Ising chain in magnetic field: 

Hmci = -JJ2f=i o-J-iO-JcrJ+i + ^Y^f=i ^J-i^J + ^Y^f=i ^^^^ '^^^ '^^^ apply GA in 

a;-basis and GGA in y- or z- basis. The conditions of the statement again hold if one properly divide 
the total Hamiltonian into the free part and the interaction. 

— an alternate XY model: 

Haxy = Jx J^iLo ^2i+i^2i+2 + ■^y X^iLi ^2i^2i+i- ^ud eq. can be applied in the z-basis. 
Particularly interesting, this model has an highly degenerate ground state, which might lead it to 
violate the third law of thermodynamics and thus to a non-zero value for ^^g^^^ \ t=o- However one 
in generals expects the zero-point entropy to depend only on the degree of degeneracy of the ground 
state, while the value predicted by GA depends on the energy of the ground state, which is clear 
from eq. . Even in this case we then expect the exact solution to be very different from GA in 
the limit of low temperature. 

5 Remarks and conclusions 

One could be wondering whether GA might be more efficiently exploited to get an insight on the 
physics of statistical models when analyzed in their non-thermal parameter space. The possibility of 
detecting signatures of quantum critical points (QCPs) in the behavior of thermodynamic functions 
is of great interest in the context of quantum phase transitions (QPTs) [TH] and in connection with 
quantum correlations between sites. It is well known that an Ising chain described by p2p indeed shows 
the presence of a QCP when h = J, leading to a transition from a ferromagnetic to a paramagnetic 
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phase. Such a transition is a second order QPT, which means it leaves signatures on the critical be- 
havior of the second derivative of free energy with respect to one of the two Hamiltonian parameters 
when T — > 0. Such a criticality is due to a non-analytical point in the ground state energy (which by 
definition is the limit of free energy as T — > 0) and is thus clearly present, as easily checked, also in 
the GA which strongly depends on Sg. Nevertheless, as shown in the previous section and as can be 
seen from ([M)) . the knowledge of Eg is a necessary condition for the application of GA to any system. 
This means that in order to calculate fa one should know in advance Eg and thus its whole set of non- 
analytical points. GA is then unable to supply further information on QPTs or to extract them more 
efficiently directly from the Hamiltonian operator. Differently from what argued in [10] then, GA can 
actually be applied to the study of QPTs (indeed at T = the free energy, as given by GA, coincides 
with the true one even if, as noticed in |10] . the two partition functions actually differ) but requires 
the same computational effort as an exact analysis and is thus not useful in simplifying any calculation. 

The prime goal of our study has been the investigation on the range of validity of the GA. The larger 
is the system, the better GA should work, therefore we mainly concentrated on the thermodynamic 
limit. We have shown that in this limit GA leads to simple expressions for thermodynamic functions 
in two important cases, namely when the temperature is high enough and when a product-state basis 
{\a)} exists for which the standard deviation of the Hamiltonian Aa does not depend on \a). Our 
results can be summarized as follows: 

— At high temperature any thermodynamic quantity as given by GA coincides with the high-temperature 
expansion of the exact function, up to terms of the order of /3^. 

— For systems with Aa independent on \a) and for systems with finitely- many blocks GA predicts 
nonzero entropy at T = 0, in contradiction to the third law of thermodynamics. 

— For systems with Aa independent on \a) GA predicts a fake thermal phase transition at an in- 
termediate temperature, in contradiction to the general result that thermal phase transitions are 
absent in ID systems with finite-range interactions [5]. 

— As discussed above in this section, GA is useless when applied to the analysis of quantum criticality 
at a fixed temperature. 

In conclusion, wc have provided convincing arguments showing that GA, generally speaking, fails at 
low and intermediate temperature, and is thus a pretty weak approximation tool. 
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